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Evolution equations of YFS and DGLAP types in leading order are considered. They are 
compared in terms of mathematical properties and solutions. In particular, it is discussed how 
the properties of evolution kernels affect solutions. Finally, comparison of solutions obtained 
numerically are presented. 



1 Introduction 

Generic QCD evolution equation covering several types of evolution reads: 

dtD{t,x) = D{t,-)®'yi{-){x), (1) 

where the change of a structure function D with a process scale t is given by its convolution 
with the evolution kernel %. I will set the boundary condition D(0, x) = 5{1 — x). 

Equations of the type ([I]) appear very often in high energy physics. For the purpose of 
this work, we will consider two simplified cases, originating from QED and QCD radiative pro- 
cesses. They were chosen due to their mathematical similarities to demonstrate how properties 
of evolution kernels influence their solutions. 

Let us consider the evolution equation with an additive convolution rule: 

dtDYFs{t,x) = / dydu 6{x-y-u)XYFs{y)DYFs{t,u) = / du %yfs{x -u)DYFs{t,u) (2) 
Jo Jo 

with the kernel given by: 

1 f^^^ du 

%YFs{u,x) = e(u-x)e(l-u) 6(u-x) . (3) 

u — X J ^ u — X 

The distribution Dyfs has been calculated by Yennie, Frautschi, and Suura J^in the context 
of exponentiation infrared singularities in QEEi. 

Another example is the evolution equation of DGLAP type^, describing evolution of gluonic 
momentum distribution in the leading order: 



dz fx\ 

dtD{t,x)= / dzdu 5{x - uz)%DGLAp{u)D{t,z) = I — Xdglap [-] D{t, z), 
Jo Jx z 



(4) 



"Precisely, we consider distributions obtained from the original one with the substitution x ^ \ — x and 
u — > 1 — w.The definition of ViYPS follows from the famous photon spectrum K(a;) = ^ — (5(s) \^ ^ 



The function D{x) is by some authors denoted as xD{x) - . The QCD sum rule requires 



[ dxD{t,x) = l. 
Jo 



(5) 



The above can be fulfilled if ^ 

/ dxD{0,x) = l, (6a) 

^0 

dt [ dx D{t,x) = 0. (6b) 
Jo 

The former condition is assured by initial condition, the latter may be obtained from ([4]): 



dt / dxD{t\x) = / dx dz du5{x — uz)%DGLAp{z)D{t,u) 
Jo Jo Jo Jo 

dzXDGLAp{z) / duD{t,u) 
Jo 



(7) 



and is fulfilled if ^ 

/ dzXDGLApiz) = 0. (8) 

Jo 

The condition ([8]) is assured by the definition %dglap = ( 



2 Comparison 

Evolution equations ([2]) and ^ are governed by similar kernels, with major difference in the 
convolution type. To compare them we will introduce regularized kernels and transform additive 
evolution into multiplicative. That form would enable us to see that the difference between them 
lies in the details of kernels definitions. 

In order to find out interrelations between both evolutions we will work using regularized 
kernels. Both kernels can be expressed as a sum of regular ("0") and singular ("5") parts 

3C(x) = %^{x)Q(l - X - e) - %\x)5{l -x-e). (9) 

Parameter e <C 1 plays the role of a constant (x-independent) infrared cutoff. 

We will write both kernels in the form they are used in evolution equations ([2]) and dl]), as 
functions of two variables. 

In the YFS case: 

%YFs{u, x) = ^e(n -x-e)-6{u-x-e)\Ti-. (10) 

U — X £ 

The DGLAP kernel %dglap is usually regularized in the (x,z) plane as follows: 

•y^DGLApi.z) = -^G(l -z-e)-5{z-l-e) C ' dz-^. (11) 

I — Z Jq L — Z 

The above notation stresses that the kernels enter their evolution equations regularized 
in different variables, according to their convolution rules. In the YFS case e is defined in 
"additive" variables (x, u) and assures x — u > 0. The DGLAP kernel is usually regularized 
"multiplicatively" , such that f < 1- 

Simple transformation of (UD leads to an evolution equation resembling the YFS additive 
convolution (12]): 



dtDit,x) 



^ '"9(1 - z - e) ^ / x\ 1 ^. I ^ / ,x\ I 

^ ^ D it,-] - - 6{z - l)ln-D (t\ 



1 - z 
e{l-x/u-e) 



ZJ z 



e 



ZJ z 



dz 



1 



U — X 



D{t, u) - 6{x - u) In -D{t, x) 



(12) 



du. 



^DGLAP in p2|) . however is regularized in (x, z) variables, which means that the cutoff is u- 
dependent. This regularization is given by 0(1 — x/u — e), which transforms into u > -j^. From 
(jl2p it is already transparent that different regularization of evolution kernels led to differences 
in their "0" parts. 

Since e is just a parameter set at the end, one can go a step further and regularize DGLAP 
kernel on the {x,u) plane with constant parameter e': u — x > e' and compare both kernels in 
that language. In that case, e in (fT2|l is x-dependent and should be expressed by x and e' . It 
could be done as follows: Let e be a (so far not specified) function of x and e'§. One can rewrite 
the last line of ()12p in this new regularization: 



dtD{t,x) 



U 



l-e(a:,£') 



U — X 



D{t, u) - 6{x - u) In — jzD{t, x 



e[x, e' 



du. 



(13) 



from which it follows that the requirement of constant cutoff e' leads to a formula 

£{x,e') 



1 — e{x,e')' 



(14) 



relating both cutoffs. From (jl4p we obtain e{x,£') = ~ — 
regularized on the (x, u) plane by constant e' has the form: 

1 X 

^DGLAp(x, u) = Q(u — X — e') 5{u — x) In — . 

u — X e' 



— so that DGLAP evolution kernel 



(15) 



Now it is clear that both kernels, although regularized identically and having equal "0" 
parts differ in their singular parts and are therefore different distributions. 



3 Normalization properties 

Another remark concerns normalization of both kernels. It is an intrinsic property of ker- 
nels, independent of regularization techniques and following from definitions. %dglap is nor- 
malized according to ([8]), whereas yivps inherits its normalization from K(?;). J^'K.{v)dv = 
¥<.{u,x)dx = 0, whereas the integration range for %yfs{u-,x) is [0,ii]. Therefore 

/ %YFs{u,x)dx= I K(w)du = alnn < 0. (16) 
Jo Jo 

As %YFs{u, x) is not normalized according to evolution equation, the normalization of Dyps 
is time-dependent (and decreasing). Normalization properties of %dglap are imposed by ([5]), 
so D DGLAP is normalized by construction. 



''Redefinition e — > e(x) is a legal operation, as x is just a parameter of the integrand (|12p . 
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Figure 1: Numerical results 

4 Comparison of solutions 

The graphs in Figure |4] present solutions of ([2]) and ^ with common initial condition D{tQ = 
0,x) = 6(1 — x), obtained at chosen evolution times. 

For small evolution times, both distributions are almost identical. The best agreement is in 
the region x ~ 1, where both (|1U|) and (|15p are very similar. In QCD D(t\x ~ 1) represents 
contribution to the total parton density function coming from very few/ very soft emissions 
that do not change the initial momentum of the branching gluon. In YFS evolution this region 
corresponds to y ~ in its physical variables meaning that emitted photons are "reasonably 
soft" and do not break the energy conservation. For small t one solution may be approximated 
by the other. 

The discrepancies arise for small x, due to differences in singular parts of evolution kernels. 
They grow larger with increasing t. 

For t — > oo, DGLAP solution converges to 6{x = 0) , and YFS has a non-singular asymptotics 
at X — > 0, ~ x"*. The right graph in Figured] presents comparison of solutions for large evolution 
time. It is visible, that the normalization of Dyps decreases, whereas the normalization of 
D DGLAP is constant. 
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